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Abstract. The influence of a layered aperiodic modulation of the couplings on 
the critical behaviour of the two— dimensional Ising model is studied in the case of 
marginal perturbations. The aperiodicity is found to induce anisotropic scaling. 
The anisotropy exponent z, given by the sum of the surface magnetization scaling 
dimensions, depends continuously on the modulation amplitude. Thus these systems 
are scale invariant but not conformally invariant at the critical point. 



The critical behaviour of quasiperiodic or aperiodic systems is better understood 
since Luck recently proposed a relevance-irrelevance criterion [1]. As in the Harris 
criterion for random systems [2], the strength of the fluctuations of the couplings, 
on a scale given by the correlation length, is of primary importance for the critical 
behaviour. Thus an aperiodic perturbation can be relevant, marginal or irrelevant, 
depending on the sign of a crossover exponent involving the correlation length 
exponent of the unperturbed system v and the wandering exponent u> which governs 
the fluctuations of the aperiodic sequence [3]. The criterion explains earlier results 
(for references see [1]) and has been confirmed in recent works on the layered 2c? Ising 
model [4-6]. 

In this Letter, we report on some recent results supporting the occurence of 
anisotropic scaling in the two-dimensional layered Ising model with a marginal 
aperiodic modulation of the exchange interactions. 

We consider a system with a constant interaction K 1 along the layers and 
apcriodically modulated interactions K 2 {k) (in k B T units) between neighbouring 
layers at k and fc + 1. In the extreme anisotropic limit, K 1 — > oo, K 2 {k) — > 0, the 
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row-to-row transfer operator involves the Hamiltonian of a quantum Ising chain [7] 

n = ~lJ2H + x ^M + i] (i) 

k 

where the cts are Pauli spin operators and the coupling A fc is given by the the ratio 
—2K 2 (k)/\n.(tanhK 1 ). For the aperiodic system, we use the parametrization \ k — \r^ k 
where f k takes the values or 1 given by an aperiodic sequence which is constructed 
through substitution. 

In the following, we shall consider: 

• the period-doubling sequence [8] with the substitutions 5(1) = 1 0, 5(0) = 1 1, so 
that, after n iterations, one obtains 

n = 1 

n = 1 10 

(2) 

n = 2 10 11 

n = 3 10 1110 10 

• the paper-folding sequence [9] with the two-digit substitutions 5(1 1) = 1 1 1, 
5(1 0) = 1 1 0, 5(0 1) = 1 1, 5(0 0) = 1 0, leading to 

n = 11 
n = 1 110 1 

(3) 

n = 2 110 110 1 

n=3 1101100111001001 

• the three-folding sequence [10] which follows from the substitutions 5(0) = 010, 
5(1) = 1 1, giving 

n = 
n = 1 10 

(4) 

n=2 010011010 

n=3 10 110 10 10 110 110 10 110 10 

Most of the properties of a sequence can be deduced from its substitution 
matrix [11] with entries given by the numbers of digits (or pairs) of the different 
types in the substitutions. 

On a chain with length L, the fluctuations of the couplings can be measured 
through the cumulated deviation from their average A, which behaves as [3] 



A(L) = 2(A fc -A)«*L"F( 
fe=i ^ 



InL 
lnA x 



(5) 



Here 5 = A(r — 1) is the amplitude of the modulation. The wandering exponent w is 
given by the ratio 
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where A x is the largest and A 2 the next-to-largest eigenvalue of the substitution 
matrix. F{x) is a fractal periodic function of its argument with period unity [3]. 
Under a change of the length scale by b=L/L', the amplitude S is changed into [4] 

5' = b UJ - 1+1 / u S , (7) 

where v is the correlation length exponent. This leads to the Luck criterion [1, 12, 
13] according to which the aperiodic perturbation is relevant (irrelevant) when u > 
(<) 1 — 1/za On the border, where iv = l — 1/v, the perturbation is marginal and leads 
to (5-dependent exponents. For the 2d Ising model with v=l, marginal behaviour is 
expected for u — which is the value of the wandering exponent for the three sequences 
mentioned above. 

The surface magnetization of the Ising quantum chain takes the simple form [14] 



-1/2 

i + EIK 2 ] X _Z+*'' t=l-(VA) 2 . (8) 

j=i fc= 



Here A c = r~ po ° is the critical coupling [15, 4] and is the asymptotic density of the 
digit 1 along the sequence. The surface magnetization exponent /3 a = x ms , where x ms 
is the scaling dimension of the surface spins, since v=\. The surface magnetization 
can generally be evaluated recursively [4] and the critical exponent is obtained using a 
finite-size-scaling method [16]. This has been done for the period-doubling sequence 
for which [4] 

In (\l /2 + Xc 1/2 ) 

x ms=~ s = C 2m2 C . K = r ~ 2/3 > (period-doubling) . (9) 

where a^TJ is the exponent on the surface corresponding to the right-hand end of the 
sequence. The period-doubling sequence in (2) is symmetric if one ignores the last 
digit which does not influence the critical behaviour. As a consequence x ms = a^T,, 
i.e. the surface exponents are the same on both sides. 

Similar calculations for the two last sequences give the following results [17] 

ln(l + X 2 r ) ln(l + A- 2 ) I/, , , , 

X ™ = 2ln2 ' ^ = 21n2 ' ^ c = ' (paper-folding) , (10) 



^ = ln( 2m3 C " 2) ' ^ J ^^> K = r~V\ (three-folding), (11) 

where the values for the right-hand surface are obtained by changing r into r _1 since, 
except for the last digit, viewed from the right side the sequences in (3) and (4) 
are obtained by exchanging and 1. Similar expressions, involving K l and K 2 , are 
obtained on the corresponding 2d classical systems [18]. 

The excitation spectrum of the Hamiltonian (1) has been studied numerically 
for the three marginal sequences given above*. A quadratic fermion Hamiltonian is 

* Although the structure of the spectra of aperiodic Ising quantum chains has already been 
studied [19-21, 1], either the aperiodicity was an irrelevant perturbation or the scaling behaviour 
was not discussed. 



Figure 1. Log-log plots of the first excitation energy ei versus the length L of the 
chain for the paper-folding sequence with periodic boundary conditions and different 
values of r. The slopes — z are taken for values of L equal to 2 n (0)> which correspond 
to a constant amplitude. 



first obtained via the Jordan- Wigner transformation [22] , which is diagonalized using 
standard methods [23]. At the critical point, the low lying fcrmion excitations e n are 
found to scale with the size of the system L as 

t n ~ L ~ Z > Z = X ms+~s- (12) 

The same behaviour is obtained for free and periodic boundary conditions [17, 18, 24]. 

The oscillations around the power laws, as shown for the paper-folding sequence on 
figure 1, are due to a periodic prefactor which, like F(x) in equation (5), is a function 
of In Lj In A x . Oscillating amplitudes are obtained for other critical quantities as well. 
When the size of the system goes to infinity, L is replaced by the correlation length 
£ ~t~ x near the critical point, and the argument of the fractal function involves the 
ratio hit/ In A r 



Table 1. Extrapolated finite-size estimates for the exponent z, obtained from the 
fcrmion excitations e n for the period-doubling sequence. 



r 




C2 


£3 


£4 


£5 


£6 


expected 


5.0 


1.19834(5) 


1.19836(1) 


1.198357(1) 


1.19836(1) 


1.19836(1) 


1.198355(5) 


1.198356 




1.19836(2) 


1.19835(1) 


1.19835(2) 


1.198357(5) 


1.198359(4) 


1.198364(4) 




4.0 


1.14885(1) 


1.14884(2) 


1.148842(4) 


1.148845(2) 


1.14884(1) 


1.148844(4) 


1.148844 




1.14884(1) 


1.14884(1) 


1.148842(3) 


1.14884(1) 


1.148844(3) 


1.14884(1) 




3.0 


1.09465(1) 


1.094648(4) 


1.094649(4) 


1.094647(3) 


1.09465(1) 


1.094651(2) 


1.094649 




1.09465(1) 


1.094654(2) 


1.094653(3) 


1.094647(3) 


1.094648(4) 


1.094647(3) 




2.0 


1.03817(5) 


1.03817(2) 


1.03820(5) 


1.003817(4) 


1.03817(4) 


1.03817(1) 


1.038170 




1.0381(1) 


1.03817(1) 


1.038172(6) 


1.038174(5) 


1.038172(3) 


1.038171(2) 




0.5 


1.0381(2) 


1.03814(5) 


1.03815(5) 


1.03815(5) 


1.03816(2) 


1.03816(1) 


1.038170 




1.0381(1) 


1.03816(3) 


1.03817(3) 


1.0381(1) 


1.0381(1) 


1.03817(1) 





Figure 2. Anisotropy exponent z as a function of the strength r of the aperiodic 
modulation for the period— doubling (full line), three-folding (dashed line) and paper- 
folding (dotted line) sequences. 



In table 1 we give finite-size estimates for z, supporting (9) and (12), which were 
obtained from sequence extrapolation using the BST algorithm (see [25]). Chains of 
size L = 2™ + 1 up to n = 20 with free boundary conditions were used. For a given r, 
the first line in table 1 refers to data with n even, while the second line refers to n odd. 
Even and odd values of n give a different amplitude since the period of the fractal 
function is 2 in this case. 

The behaviour of the excitation spectrum in (12) is typical of a strongly anisotropic 
system [26] with a correlation length exponent — zv in the time direction ( i.e. along 
the layers). In the transverse direction, the correlation length exponent keeps its 
unperturbed value v = 1 since otherwise the perturbation would not remain marginal. 
Although anisotropic critical behaviour is rather common, what is remarkable here 
is the occurence of a continuously varying anisotropy z — z(r) (see figure 2). The 
anisotropy is also implicit in Luck's finite-size calculation of an effective sound 
velocity [1]. 

The singular part of the bulk energy density is found numerically to scale like L~ z 
at the critical point so that its scaling dimension is given by 



where t is the deviation from the critical coupling as defined in equation (8). From (14), 
the specific heat exponent is given by 




(13) 



in agreement with anisotropic scaling for the bulk free energy density [26] 



f(t,L)=b-^f(b^t,L/b) 



(14) 



a 



= l-z = 



l-x. 



m s 



- X, 



ms ' 



(15) 



a relation which is indeed verified numerically for the period-doubling sequence. It 
takes a negative value since the surface magnetization exponents are always greater 
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than |, the pure system value: the aperiodic modulation of the coupling weakens 
the critical singularities. One has to notice that the proposed analytical expression 
for the specific heat exponent of the period-doubling sequence, which is obtained by 
combining equations (9) and (15), leads to a rts — A 2 / In 2 for a weak perturbation 
where A 2 =p (X) (l — p oc )(lnr) 2 in our notation. This disagrees with Luck's numerical 
results, which involved a supplementary scaling assumption [1]. 

On a semi-infinite system, using a finite-size scaling argument [17], the scaling 
dimensions x es of the surface energy density on the left surface can be related to the 
dynamical exponent z and the corresponding surface magnetization exponent. It is 
given by 

x es =z + 2x ms (16) 

which agrees with the numerical results. A similar expression is obtained on the right 
surface. 

According to anisotropic scaling [26] , the critical spin-spin correlation function on 
the left surface is expected to behave as 

G s (r ll ,t) = b- 2x -G s (r ll /b*,b 1 /»t). (17) 

At the critical point t = 0, the decay exponent is given by 2x ms / z — 2x ms / (x ms +'x^) 
instead of 2x ms for an isotropic system. Such a modified decay has been obtained on 
the 2d classical system [18] as well as for the quantum chain [17]. When the sequence is 
symmetric the decay exponent is equal to 1, i.e. it is the same as for the unperturbed 
system. Such a behaviour is indeed obtained with the period-doubling sequence. 

Marginal aperiodic perturbations of the Ising quantum chain have been shown to 
induce strongly anisotropic critical behaviour. The anisotropy exponent z, which is 
found numerically to be the sum of the surface magnetization scaling dimensions x ms 
and x^, varies continuously with the amplitude of the aperiodicity. The values of 
other bulk and surface exponents have been conjectured on the basis of numerical 
results and scaling assumptions. The bulk magnetic behaviour remains to be studied. 
Details will be given in forthcoming publications [17, 18, 24]. 
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